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Abstract 

The kinematical characteristics of fragments 
and hght particles observed in central highly 
fragmented nuclear collisions at intermediate 
energies are compared with the results of a 
model assuming that the initial momentum 
distribution of the nucleons inside the two 
partners of the reaction has no time to relax 
before the disassembly of the system. The 
rather good agreement between the results of 
the model and the experimental data suggests 
that multifragmentation is a fast process with 
a strong memory of the entrance channel. 



1 Introduction 

In central nuclear collisions at intermedi- 
ate energies, multifragmentation, defined as 
the emission in a short time scale of sev- 
eral species of atomic number larger than 
two, is an important process as compared to 
other decay mechanisms such as the forma- 
tion of heavy residues or fission (see for in- 
stance 1^). Such reactions are thus believed 
to be the ideal tool to study the transition 
from a liquid-like state (nuclei at normal den- 
sity) towards a gas-like state associated with 
the vaporisation of the system In such a 
picture, nuclear multifragmentation is asso- 
ciated with the excursion of the excited sys- 
tem at low densities inside the coexistence re- 
gion of the nuclear phase diagram . There, 
a volume instability, the so-called spinodale 
decomposition, leads the system to fragmen- 
tation. Such a scenario is supported by mi- 
croscopic transport calculations based on the 



nuclear Boltzmann equation and its stochas- 
tic extensions 0]. These models predict a 
compression-expansion phase driving the sys- 
tem at low density. 

Recent experimental analyses based on nu- 
clear calorimetry have claimed evidence for a 
liquid-gas phase transition through the study 
of various signals [^, ||, Some of these 
analyses make extensive use of the thermal 
multifragmentation statistical models |]^, |^ 
to prove the existence of thermal equilibrium 
^. In such a framework, charge distri- 
butions and fragment multiplicities are fitted 
with help of several parameters characterising 
the thermalised source. More involved stud- 
ies are based on 'back-tracing' techniques and 
they allow to obtain not only the mean values 
but also the distributions of the parameters of 
the source |]I0[. It is worth noting that, in all 



cases, not all data are fitted. Indeed, it is very 
often assumed that a strong pre-equilibrium 
phase is present and this is why kinematical 
cuts are applied to the data before compar- 
ison with the predictions of the models. In 
such thermal models, the only source of fluc- 
tuations regarding the kinetic energy distri- 
butions of the emitted species is the thermal 
motion inside the thermalised source associ- 
ated with a possible (usually assumed self- 
similar) collective motion |lT[. Recently, de- 
formation has been introduced to account for 
the observed anisotropics in the data . 



Our aim in the present work is to ex- 
plore the possibilities of a scenario in which 
the main hypothesis underlying the statisti- 
cal models is abandoned. This means that we 
wish to compare the data (with as less kine- 
matical cuts as possible) with a model taking 
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into account explicitely the entrance channel 
characteristics of the reaction as well as the 
initial correlations of the nucleons in momen- 
tum space. This is, of course, not the first 
attempt of this kind. For instance, molecular 
dynamics models have been used to under- 
stand fragmentation in such a scenario [l^ . 
Indeed, all molecular dynamics approaches, 
whether they are 'classical' or 'quantal', con- 
clude that fragmentation is a fast process for 
which there is no time for the system to build 
compression or reach equilibrium. Here, we 
develop, in this spirit, a phenomenological 
and transparent approach allowing a detailed 
comparison with the data. Since the statis- 
tical models are able to correctly reproduce 
the 'static' properties of nuclear fragmenta- 
tion, we focus our study on the 'kinematical' 
properties for both fragments and light parti- 
cles obtained in central collisions in the Fermi 
energy range. The data used for the study 
have been collected by the INDRA collabo- 
ration near the GANIL facility \T^\. We con- 
sider central collisions for Xe-|-Sn at 32 and 
50 Mev/u |15|. This contribution is organised 
as follows: the model (hereafter called the 
ELIE event generator) is described in section 
and a detailed comparison with experimen- 
tal data is developped in section ^. Conclu- 
sions and perspectives are drawn at the end 
of the paper. 



2 The ELIE event gener- 
ator 

The main idea of the model is to consider 
the extreme hypothesis for which the initial 
momentum distribution of the nucleons has 
no time to relax before the fragmentation of 
the system and for which a (small) fraction 
of the nucleons suffers hard nucleon-nucleon 
collisions leading them to escape the system 
as free particles. For the (large) remain- 
ing part of the system, only a rearrangement 
of the nucleons in the final state inside the 
various emitted products is considered. As 



ents before deca; 





Figure 1: Center-of-mass angular (up) and 
velocity (down) distributions of nucleons in- 
side the fragments before decay for Xe+Sn 
at 32 MeV/u (solid histograms), Au+Au at 
40 MeV/u (dashed line) and Xe+Sn at 50 
MeV/u (dashed histograms). 



mentioned in the introduction, we are essen- 
tially interested in the kinematical properties 
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of multifragmentation. Therefore, the defi- 
nition of the partition (that is the 'chemi- 
cal' composition in terms of complex parti- 
cles and fragments present in the final state 
of the reaction) is not determined from ab 
initio calculations but we rather use the in- 
formation given by the experimental data. In 
other words, the mean multiplicity of each 
species is put 'by hand'. However, since it 
is assumed that the fragments are excited, 
the total multiplicity of detected particles is 
the sum of two components. The first one is 
associated with the 'fast' production during 
the reaction time and the other corresponds 
to longer time scales and is associated with 
the evaporation of particles by the fragments. 
An iteration on the magnitude of these two 
components is performed until a reasonable 
agreement with the data is reached. As far 
as fragments are concerned, only the mean 
and the second moment of the multiplicity 
distribution are needed. Indeed, it appears 
that the charge distribution can be repro- 
duced rather accurately by randomly choos- 
ing sequentially the size of the fragments (see 
Fig.^). Once the partition is fixed, the ener- 
getics of the process is considered. The en- 
ergy balance writes: 

Eq = Eko + Epot + Q + Ek + E* (1) 

where Eq is the available centre-of-mass en- 
ergy, is the energy corresponding to the 
few knocked-out nucleons, Q is the mass en- 
ergy balance between the entrance channel 
and the considered partition, E* is the exci- 
tation in the fragments while Ek is the total 
kinetic energy carried away by the fragments 



and the light particles. Last, Epot is the po- 
tential energy associated with the long range 
Coulomb force acting between the charged 
species of the partition. As far as light com- 
plex particles and fragments are considered, 
their composition is obtained by randomly 
chosing nucleons in the two initial Fermi mo- 
mentum spheres separated by the relative 
momentum corresponding to the beam en- 
ergy. The kinetic energy of each species is 
thus obtained by summing the momenta of all 
nucleons belonging to the considered complex 
product. The total kinetic energy obtained 
by this procedure is such that there is little 
chance to fuUfill energy conservation. Thus, 
a Metropolis procedure is developped. An 
exchange of a nucleon (a Metropolis move) 
between two species chosen randomly is per- 
formed and the kinetic energies recalculated. 
The procedure is iterated until the value e'^^ 
at step i matches the value of Ek needed 
to fulfill Eq.[l| with a given accuracy of typi- 
cally less than a percent. Then, the partition 
is ready for propagation in order to 'burn' 
the excitation energy of the fragments and to 
transform the potential Coulomb energy into 
kinetic energy. The initial momenta are given 
by the procedure described above. The initial 
positions are obtained by randomly placing 
the species inside a volume sufficiently large 
so that nuclei and particles do not overlap: 
this gives the value of Epot- The volume in 
question is not a critical parameter since the 
key quantity to be compared with the data 
is not the separate values of Epot and Ek 
but rather their sum (which is constrained by 
Eq.|I]) (see for instance for a detailed dis- 
cussion on this point). Then, the partition is 
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propagated in a similar way as in the SIMON 



event generator . In particular, the decay 
in flight of excited species is considered in or- 
der to preserve space-time correlations. The 
generated event is then filtered through the 
INDRA software filter and stored for analysis. 
Apart from the fact that the partition is to a 
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Figure 2: Light charged particle and fragment 
multiplicities and center- of-mass kinetic en- 
ergy distributions of protons emitted in cen- 
tral Xe-hSn collisions at 50 MeV/u (see sec- 
tion Black points: experimental data. 
Filled histograms: results of the ELIE event 
generator. Data have been normalized to the 
number of events as explained in the text. 



large extent put 'by hand', there remains es- 
sentially two parameters: The amount of fast 
knocked-out nucleons and the mean excita- 
tion energy in each fragment. These two pa- 
rameters are chosen in order to get the best 
agreement with the experimental data (see 
next section). Let us discuss first the question 
of the excitation energy, E*. Fig.|l| shows the 



kinematical characteristics of the nucleons in- 
side the fragments. The angular distributions 
are more and more forward-backward peaked 
as the beam energy increases meaning that 
the fragments in this approach are more and 
more deformed and thus more and more ex- 
cited. Since, in the model, the short range 
interaction among the particles is only taken 
into account by means of the Q- values, nucle- 
ons do not interact individually and therefore 
the velocity distributions shown in Fig.0 are 
not physical. They however show little evo- 
lution and thus, it is expected that the ex- 
citation energy resulting from the deforma- 
tion increases slowly with the beam energy. 
Anticipating the results of the next section, 
we have found a good agreement with the 
data for the following values oi E*\ 3.5 and 
5. MeV/u for resp. 32 and 50 MeV/u. These 
values are compatible with those found in 
JI^ . The other parameters of the simulation, 
that is the amount of first chance knocked- 
out nucleons and the initial multiplicities of 
light particles and fragments are adjusted to 
reproduce simultaneously the observed multi- 
plicities and the center-of-mass kinetic energy 
distribution of the proton as shown in Fig.|^. 
The percentage of knocked-out nucleons in- 
creases with the beam energy as it should be 
since the phase space for in-medium nucleon- 
nucleon collisions opens gradually. In terms 
of the fraction of the total number of nucle- 
ons, one finds: 2. and 4. % for resp. 32 and 
50 MeV/u. The mean and the width of the 
fragment multiplicity distribution are 6.5, 3 
and 9, 3 resp. for 32 and 50 Mev/u. 
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Figure 3: Up, left: experimental correlation 
between 9 flow o-nd Etransi2 for Xe+Sn colli- 
sions at 50 MeV/u. Up, right: same hut for 
the simulation. Down, left: distribution of 
9 flow- Black points are experimental data, the 
filled histogram is the result of the simulation 
normalised to the largest values of 9 flow- The 
open points corresponds to the difference be- 
tween the data and the calculation. Down, 
right: same but for Etransi2- Here, the nor- 
malisation used is the same as for 9 flow 



Figure 4: 3-D plot (in log scale) of the par- 
allel velocity vs charge of selected events for 
Xe+Sn at 50 MeV/u. Up, left: experimen- 
tal data. Up, right: experimental data after 
'subtraction' of the simulation (see text for 
explanation). Down, left: experimental data 
for central events as defined in the text (in 
this case, events with flow angles larger than 
30 degrees). Down, right: same as down,left 
but for the simulation. 



3 Comparison with the 
data 

The selection of the data has been performed 
using first completeness criteria, requiring 
that at least 80 % of the total charge and 
total linear momentum (for charged parti- 



cles) be detected. This is a necessary con- 
dition to perform an event by event analysis. 
Then, for each event, a sorting is applied by 
means of two variables. First, a momentum 
tensor analysis is developped. The diagonal- 
isation of the tensor gives three eigen- values 
with which several sorting variables may be 
defined. Here, we use the so-called ^//ot« an- 
gle. This is the angle between the main axis 
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of the tensor and the beam axis. For more 



details, see . Second, a variable using the 
characteristics of the light charged particles 
is built: This is called Etransu IHl defined 
as the sum of the transverse energies of all 
detected charged particles with charge lower 
or equal to 2. Fig.^ shows the correlation 
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Figure 5: Up, left: Angular center- of-mass 
distribution for all fragments emitted in cen- 
tral Xe+Sn collisions at 50 MeV/u. Up, 
right: same hut for CM mean kinetic energy 
as a function of the emission angle. Down, 
left: same hut for the charge distrihution. 
Down, right: same hut for the CM kinetic 
energy distrihution. 

between 9 flow and Etransu for an incident en- 
ergy of 50 MeV/u. A comparison between 
the model and the experimental data shows 
that the model can reproduce the largest val- 



ues of the Ofiow distribution while it misses 
the low energy part. This is tentatively in- 
terpreted by assuming that the low values 
of the flow angle distribution are associated 
with non-central collisions for which the bi- 
nary character of the reaction remains un- 
til fragmentation. Indeed, in our model, this 
possibility is not taken into account due to 
the lack of geometry of the simulation. In 
other words, our model in its present shape 
is only able to treat reactions corresponding 
to a total overlap of the two partners and 
is thus unable to treat deep inelatic reac- 
tions. From now, we consider only events for 
which the 9 flow distrihutions of the model and 
the data coincide. This corresponds to flow 
angles larger than 30 degrees for 50 MeV/u 
and 60 degrees for 32 MeV/u. All forthcom- 
ing figures have been normalised to the to- 
tal number of events fulfilling this condition. 
Fig.^ shows the atomic number-parallel ve- 
locity plot for central Xe+Sn collisions at 50 
MeV/u. The general trend of the experimen- 
tal correlation is reproduced by the simula- 
tion as shown in the lower part of the fig- 
ure. Fragments are essentially emitted at 
mid-rapidity. However, a closer look at the 
figures shows that the model underestimates 
the kinetic energy of the fragments around 
Z=20. This point will become clearer in the 
following. A tentitative 'substraction' pro- 
cedure has been applied to the data and is 
displayed in the upper part of Fig. ^. It 
consists in subtracting from the experimental 
data the contribution predicted by the model 
and associated with central collisions with the 
normalisation discussed above. It reveals the 
presence of a rather 'pure' fast component 
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Figure 6: Same as Fig. hut Xe+Sn at 32 
MeV/u. 



associated with the remnants of a projectile- 
hke and a target-hke in mid-central colhsions. 
However, for fragments with Z lower than 
10, a so-called neck emission corresponding 
to fragments located around mid-velocity is 
also evidenced. This mechanism has already 
been studied for the same system but for a 
different selection of the data 0, A de- 
tailed analysis of neck emission with such a 
technique is in progress. We now consider 
some properties of the fragments in Fig. ^ 
and p. Due to problems with the software fil- 
ter in the CM backward direction, we show 
only kinematical characteristics of particles 
emitted in the forward part of the angular 
range in the CM frame. The mean number 
and kinetic energy of the fragments ClS db func- 



Figure 7: Angular centre- of-mass distribu- 
tions of light particles emitted in central 
xe-hSn collisions at 50 MeV/u. Up-left: pro- 
tons, up-right: deuterons, down-left: tritons, 
down-right: alphas 



tion of the CM angle are rather well repro- 
duced. The multiplicity is somehow overes- 
timated. Here, we would like to point out 
that the comparison between the model and 
the experiment is made without kinematical 
cuts: In other words, the whole event is con- 
sidered. This induces strong constraints on 
the model since any departure between the 
data and the model for some observables has 
strong feed-back effects on another set of ob- 
servables. A general 'fitting' procedure such 
clS cl back-tracing analysis should be envis- 
aged in this context. The kinetic energy dis- 
tribution is nicely reproduced until large val- 
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ues showing the important role played by the 
internal motion of the nucleons inside the two 
partners of the reaction. Same conclusions 
can be drawn concerning the results at 32 
MeV/u (Fig.^. However, here, a sizeable 
deviation in the angular distribution is ob- 
served at small angles indicating the presence 
of a contribution which is not accounted for in 
the model. We will come back to that point 
later. We now discuss the characteristics of 
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Figure 8: Same as FigQ but for the mean ki- 
netic energy as a function of the CM emission 
angle and for Xe+Sn at 32 MeV/u. 



light charged particles. We recall that most 
of them are produced by picking randomly 
nucleons from the initial two Fermi spheres. 
However, some of them are evaporated by the 
excited fragments on longer time scales. Both 



contributions lead to anisotropic angular dis- 
tributions as shown in Fig.^ We stress here 
that even if the final angular distribution of 
the fragments may be isotropic (as is the case 
if one considers only events with 9 flow larger 
than 60 degrees), those particles which have 
been evaporated will show an anisotropic an- 
gular distribution. Although the proton dis- 
tribution is underestimated at very forward 
angles, the general trend for other particles is 
well reproduced. A comparison of the mean 
kinetic energy as a function of the CM emis- 
sion angle (Fig.^) exhibits a sizeable devia- 
tion, in particular for tritons and alphas, be- 
tween the model and the data. This means 
that the model is far from perfect and some 
aspects of the emission process for light parti- 
cles is missing in our formalism. It is however 
surprising that the general trends of the an- 
gular and kinetic energy distributions of light 
particles are rather correctly reproduced by 
the model using our crude 'coalescence' algo- 
rithm. 

Let us conclude this brief comparison be- 
tween ELIE and the experimental data by 
coming back to the kinematics of the frag- 
ments. In Fig.p and [T^, the mean and the 
width of the CM kinetic energy of the frag- 
ments function of their atomic number 
are displayed. At 50 MeV/u, the second 
moment is very correctly reproduced while 
the mean kinetic energy is underestimated 
as far as the charge exceeds values of 15-20. 
This could have been anticipated looking at 
the results of Fig.^: The memory of the en- 
trance channel is stronger for the data than 
for the simulation. In other words, the model 
(by means of the Metropolis moves described 
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above) allows too much mixing between the 
nucleons of the target and of the projectile as 
far as medium and large mass fragments are 
concerned. This means that an improvment 
of the model is needed to account for 'geo- 
metrical' effects. In its present version, our 
approach can only deal with those collisions 
for which an almost complete overlap of the 
projectile and the target is achieved. This 
corresponds to very small cross-sections and 
obviously does not cover the whole impact 
parameter range associated with multifrag- 
mentation. By considering only the trans- 
verse degrees of freedom (lower part of the 
figures) , one gets rid more or less of such dif- 
ficulties and the model then gives better re- 
sults. An analysis (not shown here) of the 
data at larger impact parameters shows that 
the transverse energies of the clusters are al- 
most independent of the centrahty of the col- 
hsion except for 'controlable' final state in- 
teraction. Our model should be able to re- 
produce such trends. Work on this point is 
in progress. 

At 32 MeV/u, a large deviation between 
the data and the model is observed for the 
widths although a rather correct agreement 
is obtained for the mean values. We think 
the deviation is due to an additional mecha- 
nism in the data that is not considered in the 
model: This is the remaining binary charac- 
ter of the collision despite the already severe 
selection on the impact parameter. Deep in- 
elastic collisions lead to a strong increase of 
the fluctuations in the kinetic energy distri- 
butions of the fragments, a fact that cannot 
be accounted for in the present version of 
our model. It is worth noting that no bet- 



ter agreement is achieved when comparing 
with tranverse quantities. This suggests that 
the disagreement is not only due to geomet- 
rical effects but to collective effects such as, 
for, instance, the influence of the mean field 
on the dynamics of the reaction. 





Atomic number Z Atomic number Z 



Figure 9: Up, left: Mean kinetic energy 
in the center- of -mass as a function of the 
atomic number for fragments emitted in cen- 
tral Xe+Sn collisions at 50 MeV/u. Up, 
right: same hut for the second moment of 
the kinetic energy distribution. Down, left: 
same as up, left but for the transverse ener- 
gies. Down, right: same as up, right but for 
the transverse energies. 
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Figure 10: Same as Fig^ hut for 32 MeV/u. 

4 Conclusions and per- 
spectives 

We have studied the kinematical character- 
istics of nuclear products emitted in central 
highly fragmented nuclear collisions at in- 
termediate energies. Data have been com- 
pared with a phenomenological event gener- 
ator, ELIE, in which it is assumed that the 
initial momentum distribution of the nucle- 
ons has no time to relax before the disas- 
sembly of the system. This distribution is 
thus used to sample the kinetic energies of 
the various species produced in the reaction 
under the constraint of conservation laws. It 
is worth noting that this technique is flex- 
ible enough to test rapidly other (possibly 



thermalised) momentum distributions when 
the system breaks. As far as partitions are 
concerned, they are produced 'by hand' ex- 
cept for the charge distributions which are 
reproduced by a combinatorial random algo- 
rithm, a puzzling result which remains to be 
understood. A rather good agreement with 
the data has been obtained. The kinetic en- 
ergy and angular distributions of both frag- 
ments and light charged particles are well re- 
produced. 

In particular, in our approach, both the 
collective energy and the deformation often 
claimed as being necessary to reproduce the 
data in the framework of thermal statistical 
models are 'naturally' obtained. 

We thus believe that our model is a valu- 
able alternative to thermal statistical ap- 
proaches based on shape and/or momentum 
equilibration in the matter at low density. In 
our approach, most light particles as well as 
fragments are produced rapidly at nearly nor- 
mal density. As such, they keep a strong 
memory of the entrance channel. In par- 
ticular, fragments emerge from the reaction 
highly deformed. Then, they decay sequen- 
tially to reach their ground state by evapora- 
tion on longer time scales by a rearrangement 
and emission of nucleons and possibly other 
complex particles. The strong memory of the 
entrance channel and the importance of the 
internal motion of the nucleons inside the two 
partners of the reaction are such that, in our 
picture, the origin of kinetic energy fluctua- 
tions is to a large extent non-thermal. We 
would like to point out that this picture is 



very similar to the one proposed in |21 



Deviations between the model and the ex- 
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perimental data are expected as far as the in- 
cident energy increases. Work along this hne 
is in progress. Indeed, nuclear transparency 
is expected to decrease as the phase space 
for in-medium nucleon-nucleon collisions in- 
creases. Therefore, the two Fermi distribu- 
tions should be progressively more and more 
relaxed on shorter and shorter time scales as 
the incident energy increases. In the near fu- 
ture, we also plan to extend the model to asy- 
metric systems and to non-central collisions. 
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